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Abstract 

We show that in modified Faddeev-Jackiw formalism, first and 
second class constraints appear at each level, and the whole constraint 
structure is in exact correspondence with level by level method of Dirac 
formalism. 
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The standard method for analysis of the constrained systems is known to 
be the Dirac formalism 0. Faddeev and Jackiw 0, however, have proposed 
an alternative method (FJ formalism). Their approach is based on solving the 
constraints at each level, and using the Darboux theorem || to find a smaller 
phase space together with a number of additional coordinates. The equations 
of motion of these additional coordinates will survive for new constraints, and 
the procedure will be repeated. 

In a newer approach to FJ formalism [|[ |5|, called symplectic analysis 
or modified Faddeev- Jackiw (MFJ) formalism, instead of solving the con- 
straints, one adds their time derivatives to the Lagrangian and considers the 
corresponding Lagrange multipliers as additional coordinates. In this way, 
constraints would be introduced in the kinetic part of the Lagrangian, rather 
than the potential. There are some efforts to show the equivalence of FJ or 
MFJ method with Dirac method. In ||, using a special representation of 
constrains, it is shown that principally first and second class constraints do 
appear at a typical step of FJ formalism. However, within the Dirac for- 
malism there exists a well-established constraint structure, such that at each 
level of consistency, constraints divide to first and second class ones |§ . Then 
the consistency of the second class constraints determines some of Lagrange 
multipliers, while the consistency of first class ones leads to constraints of 
the next level. 

Such a constraint structure is not known in FJ or FJM formalism. We 
show in this paper that the constraint structure of MFJ formalism emerges 
in the same way as in Dirac formalism. However, this structure, which is 
exactly similar to that given in [|J, is somehow hidden within the Darboux 
theorem or within general statements. It should be noted that some signals 
of this structure has been recognized in @, but that belong to cases where 
only first or second class constraints are present. 

Suppose we are given the first order Lagrangian 



in which y l are coordinates not nessesarily canonical of a 2N dimensional 
phase space. The equations of motion read 



L = a l (y)y l -H{y) 



(1) 



fijij 3 = diH 



(2) 



where di means djdy 1 and 



(3) 
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The tensor is called the Presymplectic tensor. We denote it in matrix 
notation as /. If y % are chosen to be the usual canonical variables (q,p) then 
the term cii(y)y l in (Jj]) would be p^q 1 . In other words, in some general (non 
canonical) coordinates y % of the phase space, the term ai{y)y l has the same 
role as ptq 1 . Therefore, according to construction of first order Lagrangian, 
it is reasonable that / is not singular. In fact the constraints would be 
imposed thereafter to a Lagrangian with nonsingular presymplectic tensor, 
as we will see. However, if for some artificial Lagrangian, / is singular with 
null eigenvectors v l a (y), then constrains $ a (y) = v l a (y)diH(y) would result 
from equations of motion (13). In this case, one can treat with the constraints 
in the same manner as we will do in the following, and there in no serious 
difference. So, from now on we assume that the presymplectic tensor / is 
not singular. Let denote components of / _1 as f % K Then equations (§) read 

/)' = FdjH. (4) 

The antisymmetric tensor / u defines the Poisson brackets of phase space 
coordinates: 

{.'/'•. '/'} /'•'• (5) 

Assuming that 

{F(y),G(y)} = d l Fd J G{y\y>} (6) 

the Poisson bracket of two arbitrary functions F(y), G(y) can be obtained as 

{F(y),G(y)} = d t FFd j G. (7) 

Suppose we want to impose the primary constraints ®$ (f 1 = !>•••) M) to 
the system given by Lagrangian (JI|). This can be done by considering the 
Lagrangian 

L' = ai (y)f-H(y)-v^ (8) 

with V/j, as additional coordinates. Hence the presymplectic tensor / should 
be replaced by the (2iV + M) x (2iV + M) singular tensor 




Multiplying both sides of (0) with the null eigenvevtors of @ restores the 
constraints . Therefore the additional term v^&®> in @ has no new 
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result. However, consistency of the constraints during time imposes new 
constraints on system, i.e. d^/dt should vanish. Adding the term to 
L, the modified Lagrangian would be 



where 



L« = (a i -rfA ld )y i -H(y), 
A,, = d&V. 



(10) 



It should be noted that , as shown in |J, solving the constraints and consider- 
ing the equations of motion for the non canonical coordinates are equivalent 
to considering the consistency of constraints $| t °- ) . 

Considering Y = (y 1 ,^) as coordinates, the symplectic tensor /M reads 



./ 



(i) 



/ 


A 


-A 






(12) 



where the elements of (2iV x M) matrix A are defined in and A is the 
transposed matrix. The equations of motion in matrix notation then reads 



/My = dH. 



(13) 



Using operations that keep the determinant invariant, it is easy to show 



that 



det/M = det 



/ 







Af-'A 



(det/)(deti/^ 1 A). 



(14) 



Assuming that det / ^ 0, /M would be singular if C = Af 1 A is singular. 
Using (H) and (|Tl|) we have 



a 



fJ,V 



{$(o),$(o)}. 



(15) 



Therefore, the singularity disappears if C^ v is invertible (this is the case con- 
sidered in 0). In Dirac terminology this means that all primary constrains 
are second class. In this case /M defines a new bracket in the same way 
as in (I?]). It can be shown that: 



./ 



(i)' 



/- 1 - f-'AC-'Af- 1 


-]~ X AC- 1 


C- l Af- 1 


c- 1 



(16) 
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Using ([H]) one can easily see that the new bracket between two functions of 
the original phase space is the well-known Dirac bracket: 



{F,G} DB = {F,G}-{F,^}C^{^\G} 



(17) 



where C^ v is the inverse of C, 



Now suppose that C^ v = is singular. As in the framework of 

Dirac formalism || , we can assume that the constraints divide to M' first 

and M" second class constraints and $5, respectively (M' + M" = M) 
with the following property 



{*«?,*<?} «o 







(18) 



det Cyv« 7^ 



where ~ means weak equality, i.e. equality on the surface of constraints 
known already (here, Primary constraints). Defining 



A^n z 



(19) 



the tensor /W reads 





( f 


A" 


A' 




-A" 










V -A' 









(20) 



Using (|IBD and ([15|) it is easy to show that the rows of 

(i'/^Ul) (21) 

which is an M' x (2iV + M" + M') matrix, are null eigenvectors of Mul- 
tiplying both sides of ( |T3D with each row of (CT), say row gives the new 
constraint 

*$> = {*$\H} = 0. (22) 

This is exactly the same result as obtained in Dirac formalism |J, i.e. the 
consistency of first calls constraints at each level gives the constraints of the 
next level. 
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On the other hand, note that the matrix in ( p0[ ) has an invertible 
sub-block 



./i 



(i) 

inv 



/ 


A" 


-A" 






(23) 



The inverse of (|23|) is similar to fllCf ). Dividing Lagrange multipliers rf to rf' 
and rf corresponding to first and second class constraint respectively, one 
can use the inverse of ( p3|) to find ?) M from ([13|) . The result is 



if 



(24) 



where C M u is the inverse of 



{•SUiS}. 



(25) 



Adding a total derivative, one can replace the term r]^"^/, in Lagrangian 
( [TD| ) with —77^"$^/ . Then from (^), one observes that the replacement 



(26) 



suffices to omit second class constraints in the remaining and to work 

only with first calls constraints . Hence, before going to next level, we 
assume that our Lagrangian is 



(27) 



where all primary constraints are first class. (These are in fact the 

previous $S , where for simplicity we have omitted the primes.) The corre- 
sponding symplectic tensor reads 



(i) 





A (0) 


_i(o) 






(28) 



where A^ = stands for gradient of primary constraints which, as 

stated, are assumed to be first class so far. 
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Now let's go to the level step. This aim would be achieved by adding the 
consistency term — t^C 1 )/^^^ ) to the Lagrangian ( p7|) to give 



L (2) 



, * <9$ (0) , * <9$ (1) ' 



(29) 



Taking m 1 ' as new variables, we have 



./ 



(2) 



^0) 



AW \ 



(30) 



/ 



As we will see (for an arbitrary level), det vanishes if the matrix Cj^ 



[IV 



c^Wj is singular, and vise versa. If det ^ then the singularity 
disappears at all. This means that the system is completely second class. If, 
however, is singular, then one can in principle divide $W to <§>ty and 

$5, and $(°) to and in such a way that 



Rearranging one can find the invartible sub-block: 



(31) 



f- 



(2) 
inv 



/ 




AW" 















o , 



(32) 



Using inverse of (|32|) , fj^) and can be derived from the equation of motion 



QTB"!). This leads to defining a newer Hamiltonian i?( 2 ) whose explicit form is 
not important. Moreover, multiplying the equations of motion with the null 
eigenvectors of 

(IW 7 (-i) )0) 0,0,l) (33) 
gives the third level constraints 



(1)1 



(34) 
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To complete the discussion, and see some more details, let us consider 
the n-th step. 
Suppose 

$W = (35) 

are n-th level constraints; 



(n-l 
k=i 



is the n-th level Lagrangian and 

( f 

f(n) = 



j yi _ H n {y) 



A(°) ■ ■ ■ \ 



(36) 



_ A (0) 

















(37) 



is the n-th level symplectic tensor, where 

At = dM!? ] . 



in 



(38) 



Omitting second class constraints up to level (n — 1), we have assumed that 
constraints $[ t °- ) , $[ t 1 - ) • • • are first class. Using operations that leave the 

determinant invariant, one can show that 



/ P 



det/ (n) = (det/).det 



01 



V 



pOn \ 
pnn 



where 



pki = f(-i) A i = |$(fc) > 3(0} . 
Using Jacobi identity and the recursion relation 

it is easy to show 



(39) 

(40) 
(41) 

(42) 
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Since $W . . . $(« i) are assumed to be first class, using ( f4"2|) we observe 
that 

P fc/ « if k + l <n 



(43) 



In other words the matrix P has the following form 

/ 
••■ p 1 ^- 1 ) 



pOn 
pin 





pnO 



p(n-l)l 
pnl 



p(n—l)(n—l) p{n—l)n 



pn(n—l) 



pi. 



(44) 



Then using (42) we can write 

det P oc (det P 



On 



(45) 



Therefore, det / (n) vanishes if is singular. This is also an impor- 

tant observation in Dirac theory ||, that the criterion to distinguish between 
first and second class constraints is singularity of their Poisson brackets with 
primary constraints. If det P 0n ^ then the singularity of the theory disap- 
pears and all of the constraints would be second class. If, however, P 0n is sin- 
gular, one should divide all sets of constraints to first class ones: $> ®) , • • • 

and second class ones: 

$ i"> • ■ ■ The matrices A®, ■ ■ ■ A^ also break 

into A^', ■ ■ ■ A^' and A^", ■ ■ ■ A^" accordingly. Rearranging such that 
A(°>", ■ ■ ■ A^ n '" come first, we will find the invertible sub-block as 



Si 



mv 





( * 


A(°)" . . . 






Jn) _ 
J inv 


_^(0)// 

















J 




found to be 










- 1 - f^AC^Af- 1 


-f^AC- 1 




C^Af- 1 


c- 1 



(46) 



where A denotes the matrix: 

A = ( A (0) " | | • • ■ A {n) " ) . 



(47) 



(48) 
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This enables us to solve the equations of motion to find rj^", ■ ■ ■V^"- Sub- 
stituting into Lagrangin, the modified Hamiltonian if ( n+1 ) would be derived. 
Finally the null-eigenvectors of (after rearranging) can be written as 

2n-l 

(IW'V^OX^,!)- (49) 

Multiplying the equations of motion with ( f49|) finally gives the constraints of 
the next level 

$ (n+l) = {$(n) H (n)^ ( 5Q ) 

As is observed, the whole thing goes on in MFJ formalism in the same way 
as in Dirac formalism, i. e. at each level a number of second class con- 
straints , • ' ' > 3^1" ) are separated, and a set of new constraints $( ( n+1 ) 
emerge. The whole procedure will finish at level N, say, in two case: first, if 
$| 1 Ar+1 ^ ) m |$^, h\ vanish weakly and next, if the symplectic matrix is 
invertible. In the first case the system possesses gauge invariance generated 
by ($^°- ) , ■ • • , (assuming that we have put away second class constraints 

in previous levels). In the second case the system is completely second class 
and there is no gauge invariance. 

As a final point, we give a few words about a seemingly contradiction 
between Dirac and MFJ methods appeared in || when discussing the model 

L = (q 2 + q 3 ) qi + q±qz + ■= (ql - 2g 2 <? 3 - ?f) • (51) 



This model was first introduced in |]TD| as an example that the equation of 
motion can not be derived from a second order Lagrangian. Then consider- 
ing qi as coordinates of a configuration space, [ll] has discussed constraint 



structure of the model in Dirac theory. In this way four second class con- 
straints would emerge and the dynamics on the constraint surface goes on 
via Dirac brackets. The model then is analysed in || using MFJ method and 
no constraint has been encountered. 

To see what has happened, we consider the general form of a first order 
Lagrangian as given by ([I]). We emphasize that in FJ method the variables 
of the first order Lagrangian are viewed as coordinates of a phase space. It is, 
however, possible to consider those variables as coordinates of a configuration 
space. Hence, suppose we are given the Lagrangian 

L = ai{q) qi - V{q) (52) 
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which is obviously singular. Following Dirac method, the primary constraints 
are 

$i = Pi - di(q) « (53) 

We see that dublicating the space of q % to phase space of (q l ,Pi) is compen- 
sated by the same number of constraints 

Since there is no quadratic term in Lagrangian (p2[) , the canonical Hamil- 
tonian is the same as V(q) and the total Hamiltonian reads 

H T = V(q) + \ i (p i -a i (q)) (54) 

where A* are Lagrangian multipliers. It is easy to see that 

{*'•*>} =a?-a? s/ «- (55) 

Suppose / is not singular. In the context of MFJ method, this leads to 
solving equations of motion (^j) to get 

q* = pdjV. (56) 

However, in the context of Dirac theory, nonsingularity of / means that the 
constraints $i are second class and consequently the constraint surface given 
by Pi = CLi{q) has a phase space structure by itself. 
The consistency conditions $j ~ gives 

{$uH T } = -d i V + f ij \ i (57) 

which determines all A J in terms of coordinates. Inserting in fl53|) gives 

H T = V(q) + (p l -a t (q))pd,V (58) 

The canonical equations of motion for q % give 

Q l = 1 T L = PdjV; (59) 
dpi 

the same result as (|56|). (The equation p» = — give time derivative of 
(|53[) , which is obvious.) 

We see that for a regular first order Lagrangian the Dirac procedure (after 
turning around) gives the same result as one can obtain directly form MFJ 
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formalism. It is not difficult to see that the same result is correct if the system 
is singular. So the apparent contradiction between Dirac and FJ methods is 
not in fact a contradiction. The basic point is considering y % (in first order 
Lagrangian) as phase space or configuration space coordinates. 

Concluding, we think that FJ and MFJ formalism are basically the same 
as Dirac formalism formulated in the language of first order Lagrangians. 
From another point of view our work here was based on extending the sym- 
plectic tensor at each step and studying it's singularity properties. This, 
however, is closely similar to method given in ||, where for a second or 
higher order Lagrangian the Hessian matrix is extended at each step and it's 
null eigenvectors are searched for. 
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